Determining the Higgs spin and parity in the diphoton decay channel 
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We calculate the diphoton distribution in the decay of arbitrary spin-0 and spin-2 bosons produced 
from gluon fusion, taking into account the fact that gluons inside an unpolarized proton are generally 
linearly polarized. The gluon polarization brings about a difference in the transverse momentum 
distribution of positive and negative parity states. At the same time, it causes the azimuthal 
distribution of the photon pair to be non-isotropic for several spin-2 coupling hypotheses, allowing 
one to distinguish these from the isotropic scalar and pseudoscalar distributions. 
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Last year July it was announced that a new boson 
with a mass around 125-126 GeV was observed by both 
the ATLAS [1] and CMS [2] collaborations. An excess of 
events was observed in 77, ZZ* and WW* production 
from proton-proton collisions at a center of mass energy 
of 7 and 8 TeV. The observed excess is consistent, within 
uncertainties, with the production and decay of the Stan- 
dard Model (SM) Higgs boson. 

Now that the existence of a new particle has been es- 
tablished, both collaborations have begun the determina- 
tion of its spin and parity. Both ATLAS [3-5] and CMS 
[6, 7] set approximately 3 a exclusions on the J p = 0~ 
scenario using the ZZ* channel and the 2+ hypothesis 1 
starts to be disfavored at the 1-3 a level in the 77, ZZ* 
and WW* channels. As the decay of a pure spin-1 state 
to two photons is not allowed according to the Landau- 
Yang theorem [8, 9], the 77 channel is being used to 
distinguish between spin-0 and spin-2 only. In the ZZ* 
and WW* decay channels, the spin-1 option should also 
be considered. 

Even though the number of events is much larger in the 
77 channel, the ability to distinguish spin-0 from spin-2 
is not much better than in the ZZ* channel. The rea- 
son is that in the 77 channel only the distribution of 
the polar angle 6 is considered [3, 10]. The spin-0 and 
spin-2 hypotheses are not very different in this variable 
after experimental acceptance cuts [3] , leading to a small 
discriminating power. The determination of the parity 
using only this angle is even impossible, as the distribu- 
tions of + and _ are exactly equal and the same holds 
true for the 2f scenarios 2 . 



In this letter we demonstrate that one can also differen- 
tiate between the different spin scenarios in the 77 chan- 
nel, by studying the dependence on the azimuthal angle 
</> in the Collins-Sopcr frame [13], which is the diphoton 
restframe with the iz-plane spanned by the 3-momenta 
of the colliding protons and the x-axis set by their bisec- 
tor. Moreover, different spin-2 coupling hypotheses that 
have an equal dependence can be distinguished from 
each other using the <j> distribution, enhancing the an- 
alyzing potential of this channel. Apart from that, wc 
update predictions for the transverse momentum distri- 
bution [14, 15] which can, in principle, be used to dis- 
tinguish the different parity states + from 0~ and 2^~ 
from 2^ in the 77 channel. Azimuthal angular distribu- 
tions have been discussed for spin-0 and spin-2 "Higgs" 
production from vector-boson fusion [16-18], but not yet 
from gluon fusion and not including linear polarization. 

A non-trivial 4> distribution in the decay of spin-2 
bosons produced from gluon fusion can be caused by the 
fact that gluons in an unpolarized proton are generally 
linearly polarized. The degree of gluon polarization can 
be calculated using perturbative QCD (pQCD) for trans- 
verse momentum of the gluon much larger than the pro- 
ton mass and is found to be large. For small transverse 
momentum pQCD cannot be used to calculate the degree 
of polarization, but this lack of knowledge turns out to 
be of little influence on the final (f> distribution, which is 
mostly dominated by the perturbative part. 

The effects of gluon polarization can be described 
in the framework of Transverse Momentum Dependent 
(TMD) factorization. In that framework, the full pp — > 
77 A cross section is split into a partonic gg — > 77 cross 
section and two TMD gluon correlators, which describe 



1 The coupling of a spin-2 boson to gauge bosons can be realized 
in multiple ways. We will use the standard notation in which 
2„ denotes a spin-2 boson with minimal (lowest dimensional) 
coupling, which is uniquely defined. 

2 The h subscript indicates that the spin-2 boson couples through 



a higher-dimensional coupling. There arc multiple highcr- 
dimcnsional couplings possible. We follow the convention of Rcf. 
[11, 12] for 2*. 



the distribution of gluons inside a proton as a function of 
not only its momentum along the direction of the proton, 
but also transverse to it. More specifically, the differen- 
tial cross section for the inclusive production of a photon 
pair from gluon-gluon fusion is written as [19, 20], 
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with the longitudinal momentum fractions X\ = 
q ■ P2/P1 ■ Pi and x 2 = q • Pi/ Pi • P 2l Q the momentum 
of the photon pair, M. the gg — > 77 partonic hard scat- 
tering matrix element and $ the following unpolarized 
proton gluon TMD correlator, 
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pi^pv jp.pi ] w hcrc P and P' are the momenta of the 
colliding protons and M p their mass. The gauge link 
U7L X in the matrix element runs from to £ via minus 
infinity along the direction n, which is a time-like dimen- 
sionlcss four-vector with no transverse components such 
that C 2 = (2n-P) 2 /n 2 . In principle, Eqs. (1) and (2) also 
contain soft factors, but with the appropriate choice of 
£ (of around 1.5 times the hadronic center of mass en- 
ergy), one can neglect their contribution, at least up to 
next-to-leading order [20, 21]. The renormalization scale 
should be chosen around the characteristic scale of the 
hard interaction. The last line of Eq. (2) contains the pa- 
rameterization of the TMD correlator in terms of the un- 
polarized gluon distribution ff (x,p r , £, A 1 ), the linearly 
polarized gluon distribution h 1 9 (x,p T ,£,/i) an d Higher 
Twist (HT) terms, which only give 0(1/Q) suppressed 
contributions to the cross section, where Q = \/q 2 ■ 

The general structure of the differential cross section 
for the process pp — > 77X is given by [22] 
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+ Fi(Q, 9) C \wi hi 9 hi 9 \ cos(#) + O I ^ I , (3) 

where the Fi factors consist of specific combinations of 
gg — > Xq,2 — > 77 hclicity amplitudes, with F$^ involving 



amplitudes with opposite gluon helicities. The convolu- 
tion C is defined as 

C[wfg] = jd 2 p T fd 2 k T 5 2 (p T + k T -q T ) 

w(p T ,k T )f(xi,p 2 T )g(x 2 ,k 2 T ) (4) 
and the weights appearing in the convolutions as 
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The TMD distribution functions contain both per- 
turbative and non-perturbative information. The tails 
(p T >• M p ) of the distribution functions can be calcu- 
lated using pQCD, but the low p T region will inevitably 
contain non-perturbative hadronic information. To get a 
description over the full p T range one needs to extract 
the TMD distribution functions from experimental data 
[22, 23]. 

To make numerical predictions we will use a functional 
form for the unpolarized gluon TMD which has, in ac- 
cordance with the pQCD calculation, a 1/p 2 tail at large 
p T and resembles a Gaussian for small p T1 
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Preferably one would fit the parameters in Eq. (6) to ac- 
tual data, but since those are currently not available we 
will instead fit to the Standard Model Higgs boson trans- 
verse momentum distribution obtained by interfacing the 
POWHEG [24-26] NLO gluon fusion calculation [27] to 
Pythia 8.170 [28, 29], assuming a Higgs mass of 125 GcV 
and a collider center of mass energy of 8 TeV. Pythia 
does not take into account effects of gluon polarization, 
so we fit the data by setting the linearly polarized gluon 
distribution equal to zero. In this way the TMD predic- 
tion without gluon polarization agrees with the Pythia 
prediction. We think this is the most realistic choice we 
can make, because Pythia is tuned to reproduce collider 
data well. Our Gaussian-with-tail Ansatz is able to ad- 
equately fit the Pythia data, as is shown in Figure 1. 
The fit results in the following values for the parameters 
a = 38.9 GcV, a = 0.555 and M = 3.90 GeV. We are 
not concerned about the overall normalization, as we will 
be only interested in distributions and not the absolute 
size of the cross section. 

The linearly polarized gluon distribution will be ex- 
pressed in terms of the unpolarized gluon distribution 
and the degree of polarization V, i.e., 
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FIG. 1. Plot of q T C[flf[] (line) and the Pythia Higgs da/dq T 
distribution for M h = 125 GeV at s/s = 8 TeV (points). 
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FIG. 2. Plot of the degrees of polarization V m i n , V and P ra „ 
at x = Mh/y/s, with M h = 125 GeV and y/s = 8 TeV. 



such that \V\ = 1 corresponds to h x 9 saturating its up- 
per bound [30] and with the correct power law tail as first 
calculated in [19]. Calculations of the gluon TMD distri- 
butions using the Color Glass Condensate model predict 
maximal gluon polarization for large p T and small x [31]. 
Ideally one extracts the degree of polarization from data, 
but this is currently unfeasible. 

Perturbative QCD can be used to calculate the large 
p T tails of the TMD distributions in terms of the collinear 
parton distribution functions as has been done in Ref. 
[21] for the unpolarized distribution and Ref. [19] for 
the linearly polarized gluon distribution. We will follow 
a similar approach, but keep finite £ instead of taking 
the £ — > oo limit and calculate the degree of polarization 
to leading order in a s from the MSTW 2008 collinear 
parton distributions [32] evaluated at a scale of /j, = 2 
GeV. 

The pQCD calculation is only valid in the limit p T ^> 
M p . To model the lack of knowledge at low p T , we will 
define three different degrees of polarization Vmi n , V and 
V m ax i of which the first approaches zero at low p T , the 
second follows the pQCD prediction and the last reaches 
up to one at low p T . Other sources of uncertainty are the 
choices of the scales C, and \i and the omission of higher- 
order terms. We estimate this additional uncertainty, by 
varying the different scales, to be maximally 10% and 
model it by letting V m ax,min approach the pQCD calcu- 
lation ±10% for large p T . More specifically, we define 
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where "PpQCD is the pQCD degree of polarization cal- 
culated at £ = l.5\/s and we take p — 5 GeV. The 
resulting V m in, V and V ma x are plotted in Figure 2. 

We will consider the partonic process gg — > X 0: 2 — > 77 
where X is either a spin-0 or spin-2 boson, with com- 
pletely general couplings. For the interaction vertex we 
will follow the conventions of Refs. [11] and [12], where 



the vertex coupling a spin-0 boson to massless gauge 
bosons is parameterized as 



V[X -> V"(gi)V"(«fc)] = a iq 2 g^ + a 3 e^^, 
and for a spin-2 boson as 



(9) 
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+ (c 2 q 2 g^ + cr e^^)^, (10) 



where q = qi+q 2 and q = q\ — q%. The coupling to gluons 
can be different from the coupling to photons, but to keep 
expressions compact we will consider them equal. 

For the gg — > Xq — > 77 subprocess, the non-zero F 
factors in Eq. (3) read 

i^i = 16| ai | 4 + 8|a 1 | 2 |a 3 | 2 + |a 3 | 2 , 

F 2 = 16| ai | 4 -|a 3 | 4 , (11) 

and for the gg — > X 2 —¥ 77 process one has 

F^lSA+lc^st + A+'il-Scl) 2 

+ || Cl | 4 (28c 2e + c 4e +35), 

F 2 = 9A-\ Cl \ 2 s 4 g + A-A+ (1 - 3c^) 2 , 

F 3 = 3s 2 g B" [3| Cl | 2 (c 2e + 3) + A+(3c 2g + 1)] , 

F^ = 6*iRe(cicS) [3| Cl | 2 (c 2e + 3) + A+(3c 2e + 1)] , 

F 4 = 9 S 4 |c 1 | 2 [2B+ + 4|c 5 | 2 ], (12) 

where we have defined A ± = \c\ + 4c 2 | 2 ± 4|c 5 | 2 , B ± = 
|c 1 +2c 2 | 2 ±4|c 2 | 2 , c n g = cos(n6) and sg = sin(6*). Overall 
factors have been dropped, because as said we will be 
only interested in distributions and not the absolute size 
of the cross section. Unlike the case for Higgs production 
from linearly polarized photons [33], there is no direct 
observable signalling CP violation in the spin-0 case. For 
the spin-2 case there is such a clear signature, being a 
sin 2(f> dependence of the cross section, which can only 
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TABLE I. Different spin, parity and coupling scenarios. 



be present if both c\ and C5 are non-zero, implying a 
CP-violating interaction. 

In Ref. [12] a set of different spin, parity and cou- 
pling scenarios is defined. To those scenarios we will 
add 2^, and 2^„ , which will serve as examples of higher- 
dimensional spin-2 coupling hypotheses that are indistin- 
guishable in the 9 distribution, but do have a different 
distribution. The scenarios are summarized in Table I. 

In Figure 3 we show the diphoton cos 9 distribution for 
the various scenarios. Looking only at this distribution 
+ and _ are indistinguishable, as are 2^ and 2^, and 
also 2^", and 2^„. 

fdjl dq\ do- /Jd<f> dq 2 T rfcosfl do 

\1 



— 2 + 

.... o* 



FIG. 3. Plot of the cos# distribution for the various scenarios. 

In Figure 4 we show the diphoton transverse momen- 
tum distribution for the different coupling hypotheses at 
fixed 9 — 7r/2 and at zero rapidity. The positive parity 
states show an enhancement at low q T (< 15 GeV) with 
respect to the negative parity states. At high q T (> 15 
GeV) this is reversed, but with such a strongly reduced 
magnitude that it is invisible in the plot. The q T distri- 
bution can thus, in principle, be used to determine the 
parity of the newly found boson [14, 15]. Although the 
difference is small and most likely difficult to measure ex- 
perimentally, this is the only way we know to determine 
the parity in the gg — ¥ X$^ — ► 77 channel. 

Figure 5 shows the diphoton (j) distribution for the se- 
lected scenarios at fixed 9 = n/2 and at zero rapidity. 
The scalar, pseudoscalar and 2^ hypotheses show a uni- 
form (j) distribution, whereas the 2+ has a characteristic 
cos(40) dependence with an amplitude of 5.4l^g%. The 
2\, and 2\,, scenarios exhibit a weak cos(4</>) modula- 
tion with an amplitude of 1.2+°;4% and a strong cos(20) 
modulation with an amplitude of 24 ± 3% and opposite 
sign. The <j) distribution thus offers a way to distinguish 
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FIG. 4. Plot of the q T distribution for the various coupling 
schemes at 9 — it/2 and zero rapidity, using an upper limit on 
the q T integration in the denominator of Mh/2. The shaded 
area is due to the uncertainty in the degree of polarization. 



± , 2+, 2^, and 2^„ from each other, something that is 
impossible with the cos 9 distribution alone. 

We want to stress again that a sin 20 dependence im- 
plies a CP-violation coupling, which is thus very inter- 
esting to search for. Note however that Higgs bosons 
produced with positive and negative rapidity have to be 
treated separately, because those regions will have an op- 
posite sign sin 20 modulation and would otherwise can- 
cel. We also want to mention that gg — > 77 continuum 
production has a non-isotropic <f> dependence, with an 
amplitude approximately a factor 3 smaller than reso- 
nance production [22, 34], which should not be mistaken 
for a spin-2 Higgs. 

Jdql do /Jdt/j dq\ do 
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- 0^,2* 
2t 




FIG. 5. Plot of the <j> distribution for the different benchmark 
scenarios at 6 — tv /2 and zero rapidity, using an upper limit 
on the q T integration of Mh/2. The shaded area is due to the 
uncertainty in the degree of polarization. 

In conclusion, we have calculated the diphoton distri- 
bution in the decay of arbitrary spin-0 and spin-2 bosons 
produced from gluon fusion, taking into account the fact 
that gluons inside an unpolarized proton are generally 
linearly polarized. The gluon polarization brings about 
a difference in the transverse momentum distribution of 
positive and negative parity states. At the same time, it 
causes the azimuthal CS angle 4> distribution to be non- 
isotropic for various spin-2 coupling hypotheses. These 
distributions allow spin and parity scenarios to be dis- 



tinguished that cannot be done with the polar angle 9 
dependence alone. We think that these observables could 
therefore form a valuable addition to the analysis meth- 
ods to determine the spin, parity and coupling of the 
newly found boson at the LHC. 
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